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Abstract 

We investigate the long-term web structure emerging in evolutionary food web mod- 
els when different types of functional responses are used. We find that large and complex 
webs with several trophic layers arise only if the population dynamics is such that it 
allows predators to focus on their best prey species. This can be achieved using modi- 
fied Lotka-Volterra or Holling/Beddington functional responses with effective couplings 
that depend on the predator's efficiency at exploiting the prey, or a ratio-dependent 
functional response with adaptive foraging. In contrast, if standard Lotka-Volterra or 
Holling/Beddington functional responses are used, long-term evolution generates webs 
with almost all species being basal, and with additionally many links between these 
species. Interestingly, in all cases studied, a large proportion of weak links result natu- 
rally from the evolution of the food webs. 
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1 Introduction 



The debate over what stabihses complex food webs has been an active one for over a quarter 
of a century (McCann, 2000). While it is found that the stabihty of a randomly hnked food 
web model tends to decrease with the proportion of links and the number of species (May, 
1972), real food webs display a high degree of stability, in spite of being very complex. Clearly, 
real food webs have features that are ignored in randomly linked models. For example, the 
distribution and strength of links in real food webs are far from random. More realistic 
approaches use hnk values taken from real food webs (Yodzis, 1981), or assemble a food 
web by repeated addition of species from a large species pool that contains different types 
of species such as "plants" , "herbivores" , "carnivores" , and "top predators" (Morton & Law, 
1997; Law, 1999). Using Lotka-Voltcrra type dynamics, assembly models of this type lead to 
large complex webs. These Lotka-Volterra type models lack nonlinear effects such as predator 
saturation, which many consider important for food web dynamics. The species pool is 
usually composed of a large number of species which are generated in an ad hoc manner, 
by pre-defining the possible types and value ranges of links. A complementary approach 
to food web stability focuses on more realistic dynamics (that include predator saturation, 
for instance) of small systems consisting only of a few species. These small systems can be 
viewed as a small part of a food web. It has been found that the presence of many weak 
links and the use of nonlinear growth rates, that prevent predators from feeding successfully 
on many prey at the same time, stabilise the dynamics of these small systems by reducing 
population fluctuations (McCann et al, 1998). Again, the possible types and values of links 
are pre-deflned in these studies. 

In real food webs, the link strengths and the linkage pattern can change with time through 
invasion and replacement with related species, and through adaptive evolutionary processes 
(Thompson, 1998). It is therefore desirable to investigate food web models that implement 
such long-term changes due to modiflcation of existing link strengths and linkage patterns. 
In the following, we call this long-term change of the food web "evolution" , with the meaning 
of "development of the food web structure in time", while the "true" biological evolution is 
of course happening on the level of individuals. Our approach is different from the species as- 
sembly models insofar that the new species are not picked irrespective of the web composition, 
but are similar to existing species. 

Considering the long-term change of food webs adds a new dimension to the complexity- 
stability debate: even though species die out and are replaced with others during the course 
of time, ecosystems as a whole persist in time. Except for rare catastrophic events due to 
external causes, food webs preserve their general structure. This means that typical structural 
features of food webs, such as the total number of species or the number of trophic layers, do 
not collapse to zero, but rather fluctuate around some mean value far away from zero. This 
definition of stability is perhaps more relevant for understanding the complexity of ecosystems, 
than some of the other definitions found in the literature. Persistence in time of a complex 
ecosystem of course also requires a certain kind of stability of the population dynamics: while 
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population dynamics does not need to reach fixed points, oscillations must be small enough 
such that the population sizes of the majority of species remain positive after the introduction 
of a mutant or invading individual. If this were not so, each change in the web composition 
would lead to a collapse of the web, and complex ecosystems could never arise. 

In this paper, we therefore study evolutionary food web models and investigate under 
which conditions these models show a persistent complex structure with several trophic layers. 
The evolutionary model that we use was studied previously with specific types of population 
dynamics (Caldarelli et al, 1998; Drossel et ai, 2001). Here, we explore the effect of changing 
the functional response used in the model. We chose a representative selection of functional 
responses found in the literature and modified them such that they are suitable for many- 
species communities. Our aim here is to show that biologically realistic requirements, in 
particular the capability of predators to adjust their feeding rates and focus on their best 
prey, are essential to the construction of large stable complex webs. 



2 The model 

The governing equation for the population Ni{t) of species i is taken to have the form 

^ = Xj:N,{t)g,,{t)-j:N,{t)gj,{t) 

j j 

-diNi{t)-Y^aijNi{t)Nj{t), (1) 
j 

which naturally accounts for the four processes which lead to a change in the population num- 
ber Ni{t). The function gij{t) is the functional response, i.e. the rate at which an individual 
of species i feeds on species j; it depends on the population sizes, and its analytical form 
will be specified below. The constant A is the ecological efficiency at which consumed prey 
are converted into predator offspring. The first term describes population growth due to food 
consumption, while the second term describes population decline due to predation. The death 
rate, dj, will be assumed to be equal to 1 for all species. The last term describes direct in- 
terference competition between two predators i and j. It is needed to represent intra-specific 
interference competition, which ensures bounded solutions when the functional response is 
independent of predator density. This term also allows the incorporation of inter-specific 
interference competition into such models, which facilitates comparison with the predator 
dependent functional responses. 

On a larger time scale, species undergo changes, either by replacement with invading 
similar species, or by evolutionary change. In order to implement these changes, each species 
was characterised by a set of 10 out of 500 possible features (Caldarelli et ai, 1998). This 
representation gives a measure of similarity between species (the number of features they have 
in common) and allows for "mutations" or "invasions" by randomly replacing one feature 
of one individual with another. Scores between features were assigned in a random and 
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asymmetric fashion. They are a measure of how useful a feature is for its carrier at feeding 
on a species carrying another feature. The interactions, or scores, aij, between species were 
obtained by adding the mutual scores between all pairs of features carried by the two species. 
Positive scores indicate that the first species can feed on the second species and negative 
scores mean that the first species is consumed by the second. In the latter case the score was 
set to zero, since this information was already contained in the (positive) score where the two 
species are interchanged. Therefore all scores were such that aij > 0. The external resources 
were represented as an additional species of fixed (and large) population size, which does not 
feed on any species. The evolutionary dynamics consisted of the following steps: starting 
from some initial species configuration (usually the external resources plus one basal species 
feeding on it), one individual was picked at random, and one of its features was changed. 
Under the dynamics ((T)), this new species either died out, or added to the system, or drove 
one or several other species to extinction. When the dynamics had reached an equilibrium, 
the next "mutation" occurred, and the process was repeated. 

In the following, we list the different functional responses used in our computer simulations. 
Lotka-Volterra dynamics are recovered by taking 

gij = aijNj . (2) 

A Holling type II functional response, which implies saturation of consumption rates at 
high prey abundance, is given by 

^ 1 + i Ln, ' 

where the sum in the denominator is taken over all prey k of species i. 

More complicated functional responses for multi-species systems can only be found in the 
recent literature. Arditi and Michalski (1996) suggest the following generalised Beddington 
form: 

a-= ^'^^^ r4l 

' 1 + Efc hkNk + E/ CiiNi ' ^ ' 

where the first sum is again taken over all prey k of species z, and the second sum is taken 
over all those predator species / that share a prey with i. We chose the cu such that individ- 
uals belonging to the same species competed more strongly with each other than individuals 
belonging to different species. 

Arditi and Michalski (1996) also suggest the following ratio-dependent functional response, 
which implements the idea that predators share the prey: 

jr{i) 



with the self-consistent conditions 
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Here is the efficiency of predator i at consuming species j, hij is the relative preference of 
predator i for prey j, R{i) are the prey species for predator i, C{i) are the species predating 
on prey i, A'J^*^ is the part of species j that is currently being accessed as resource by species 

i and N^^^"* is the part of species k that is currently acting as consumer of species j. 

Finally, in a recent paper (Drossel et ai, 2001), we used the ratio-dependent expression 

^''^ ^ 6iV,(t) + Ek «..a,,/,,(t)iV,(t) ' ^ ^ 

where fij is the fraction of its effort (or available searching time) that species i puts into 
preying on species j. These efforts are determined self-consistently from the condition 

hit) = (7) 

This condition is such that no individual can increase its energy intake by putting more effort 
into a different prey. The sum in (jHl) is over all species k which are predators of j. The 
competition strength is set equal to one only for i = k, and is smaller than 1 otherwise. 
For this model, or for any other model with a functional response which depends on predator 
density, the last term in Eq.((l]) should be dropped. 



3 Simulation results: food web structure 

Let us first present results for the simplest type of population dynamics, which is Lotka- 
Volterra dynamics, Eq. Q. Starting with one species and the external resources, input at 
a rate R, we allowed for enough time for the system to evolve towards its characteristic 
stationary structure. A typical food web resulting after sufficiently long time is shown in 
Fig. All species feed on the external resources, and additionally on other species. Only 
for short transient periods of time does the web have species on the second trophic level, 
and virtually never on the third. We also performed simulations that started from a complex 
web with several trophic layers, which was stable under the population dynamics. However, 
the evolutionary dynamics caused this structure to collapse, and the resulting structure was 
again similar to that shown in Fig. ^ These results were obtained without inter-specific direct 
competition, i.e. with aij = c6ij where 6ij is the Kronecker delta defined by 6ij = for i j 
and 6ii = 1. They are still true when direct inter-specific competition is included, as in a 
different evolutionary model described by Lassig et al (2001) and also in our model if we set 
Oiij = cpijil + qij)l'2, in Eq. ((H). Here pij is equal to one if i and j share at least one prey and 
is zero otherwise, and qij is the fraction of discrete features shared by the two species. 

Our interpretation of these results is as follows. When a modified species that can feed on 
more prey than a parent species is introduced into the model, the modified species takes over. 
This is due to the linear growth rates of Lotka-Volterra systems that allow a species to feed on 
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Figure 1: Example of a food web resulting from an evolutionary model with Lotka-Volterra 
population dynamics. The arrow direction indicates the flow of resources, and the arrow 
thickness is a measure of link strength. Links are only drawn if a species obtains more than 1 
percent of its food through that link. The radius of the circles increase logarithmically with 
population size. 



each of several prey just as easily as it would feed on one prey, if only this prey was present. 
If allowed to evolve, an artificially constructed Lotka-Volterra food web will therefore tend to 
a simple web structure where each species feeds on many other species and in particular on 
the most abundant resource, the external environment. Clearly, this feature of the model is 
unrealistic as there are evolutionary trade-offs that do not allow species to feed very easily on 
an arbitrarily large number of species. In particular, species cannot normally feed on external 
resources and on many other types of prey. 

In the light of these arguments, a Holling or Beddington type functional response is more 
realistic, since it naturally limits the total amount of prey taken by a predator. Indeed, 
the need for nonlinear growth rates in the context of multispecies communities has often 
been discussed in the literature (Pimm, 1991; Rosenzweig, 1995; Vandemeer et al, 2002). 
Functional responses of the form (jH} were used by McCann et al (1998) in their investigation 
of the effects of weak links. They have the biologically meaningful feature that predators 
cannot maintain a high feeding efficiency on many prey at the same time, and they have a 
stabilising effect on the dynamical equations, at least for subsystems consisting of up to four 
species (McCann et al, 1998). Additional direct inter-specific competition can be introduced 
in a natural way by using functional responses of the Beddington form However, when 
evolutionary dynamics by modification of existing species are added to either of these models, 
we find again that complex webs with several trophic layers are not stable and that after some 
time the webs consist merely of basal species, with a structure very similar to that shown in 
Fig.m 

It appears that even for the Holling or Beddington type functional response too many 
species can feed on the same food source, and that the model does not have a mechanism 
that prevents possible links (i.e. a^j > 0) from being realised. In the models we are considering 
here, there is always a possible predator-prey interaction between any two species i and j, 
either Ojj > and aji = or aji > and Ojj = 0. Thus the ecosystems are potentially highly 
connected with one half of all inter-specific links realised. However, species usually focus 
on those prey to which they are best adapted. We therefore must introduce a mechanism 
which would limit the number of species a predator actually feeds on, with potential links 
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species 


Trophic level 


1 


2 


3 


4 


5 


6 


Number of species 


Frequency of occupation 


Lotka-Volterra without 
adjusted interactions 


40.6(2.5) 


2.3(0.1) 


40.5 


1.6 










1.000 


0.061 










Lotka-Volterra with 
adjusted interactions 


69.8(20.4) 


1.4(0.1) 


20.9 


22.1 


19.1 


7.5 


1.4 


1.0 


1.000 


1.000 


1.000 


0.987 


0.140 


0.003 



Table 1: Food web statistics for the two Lotka-Volterra models: without and with adjusted 
interactions. The results are averaged over ten different simulations (lasting 200000 iterations) 
and over the last 20000 iterations of each simulation. The quantities in brackets give standard 
deviations over the ten runs for the number of species and links per species. Only links between 
non-environment species that constituted greater than 1% of the predator's diet were counted 
in the calculation of the links per species. 

becoming active only when the preferred prey is scarce or unavailable, or after a change in 
the composition of the web. One way to (artificially) realise this feature is to allow only the 
best predators to feed on a given prey. We did this by replacing the interactions Ojj with 
adjusted interactions a-^- = a™'^''(l — (a™'^'' — aij)/5) (Caldarelh et al, 1998), with S being 
a small parameter and aj^^^ being the largest interaction against j. Negative a^^- are set to 
zero. Fig. |21 shows a food web obtained from Lotka-Volterra dynamics with this artificial 
constraint. Table 1 shows the mean number of species, of links per species, and the mean 
occupation numbers of the trophic levels for this model, for R = 1 x 10^, A = 0.1, c = 1.0 
and S = 0.2. For comparison, the results for the original Lotka-Volterra model with the 
parameters for R = 1 x 10^, A = 0.1, and c = 3.0 are also shown. Similar web structures were 
obtained including direct inter-specific competition or using Holling and Beddington forms 
with the same type of adjusted interactions. 

This shows that the ability of predators to concentrate on the prey that they are best 
suited to exploit, rather than on all possible prey, is essential for the production of realistic 
food web structures. 

Finally, let us discuss evolutionary food web models with ratio-dependent functional re- 
sponses. Ratio-dependent functional responses naturally limit the number of actual prey of a 
predator by dividing the available prey among its predators. One such ratio-dependent form, 
given by (0), was suggested by Arditi and Michalski (1996). However, it leads to such strong 
competition that only one predator species can feed on a given prey species and thus cannot 
give rise to complex food web structures. (One can conclude this directly from the published 
webs (Arditi & Michalski, 1996) or from an analytical calculation with two predator species 
and one prey species). In reality, however, different species exploit a food source in different 
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Figure 2: Example of a food web resulting from an evolutionary model with Lotka-Volterra 
population dynamics and adjusted interactions [6 = 0.2). The same conventions apply as for 
Fig. ^ except that vertical position is now proportional to the average path length from the 
species to the environment weighted by diet proportions. Notice that the population sizes 
vary much less across trophic levels than when a ratio-dependent functional response is used 
— as in Fig. El for instance. 



ways and therefore competition between individuals of different species is less than between 
individuals of the same species. 

This feature is implemented in the functional response (jHl), which was introduced by some 
of us (Drossel et ai, 2001). Furthermore, equations (jHl) and ((Zj) have several other biologically 
meaningful features: when prey is very abundant, the consumption rates saturate due to the 
term containing b in the denominator. Additionally, predators divide their effort (or time) 
among the available prey in order to obtain the maximum possible amount of food. The 
condition ((Zj) is an evolutionarily stable strategy and allows species to swiftly adjust to a 
modified situation, thus endowing the food web model with the flexibility found in real webs. 
Finally, equations (jHI) and ((Zj) are invariant if identical species are aggregated into a single 
species. This feature is shared by the other models used in this paper if the last term in 
equation ((T)) is either set to zero or present for any pair of species but not if, for instance, 
only inter-specific competition is included. 

Using the functional response ((Hj) in the evolutionary model, we rapidly obtain large and 
complex food webs with a stable structure, although there is an ongoing turnover of species in 
the system. The dynamics (0) converge quickly to a fixed point. The food webs generated by 
this model share many features with real food webs, such as the fraction of top, intermediate 
and basal species, and the mean number of links per species. This is discussed in detail 
elsewhere (Drossel et ai, 2001; Quince et ai, 2004a). 

The webs generated by this model often preserve their structure under changes in the 
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Figure 3: A model web generated using the functional response © before (left) and after 
(right) deletion of basal species 5. The two shaded species feeding on it managed to survive 
by feeding at a higher rate on other species. 

species composition. Fig. El shows a web before and after a basal species was deleted. One 
can see that in this case all other species manage to survive, with the strength of several 
links being modified, illustrating the flexibility and resilience of the model food webs. Of 
course, the deletion of a randomly chosen species is not always buffered in this perfect way. 
A useful measure of the ecosystem stability to deletion is the fraction of species that can be 
removed without causing further extinctions (Pimm, 1979). This statistic depends on the 
model parameters; for the values used to generate Fig. El it is 62%. For the deletions which 
do cause further extinctions the distribution of event sizes decays roughly exponentially with 
a characteristic size of just a few species (Quince et ai, 2004b). 

4 Simulation results: link strength distributions 

There has been an increasing realisation that food webs have a large proportion of weak links 
(Paine, 1992; Tavares-Cromar & Williams, 1996; Berlow et al, 1999; Neutel et a/., 2002), 
and that weak links tend to stabilise population dynamics (McCann et ai, 1998). For this 
reason, we evaluated the link strength distribution in our model webs. We found that for 
all functional responses that generated large ecosystems either the link strength distributions 
were skewed towards zero or a large fraction of links were zero. This is a highly non-trivial 
result as in contrast to other work on the topic of weak links, the link strength distributions 
in our model are an emergent property of the system and not put in by hand. It is a strong 
indication that weak links are the natural outcome of long-term ecosystem evolution coupled 
to population dynamics. 

In the literature there are different definitions of link strength, ranging from the biomass 
flowing through a link per unit time to the response of one population size to a small change 
in another population size (Laska & Wootton, 1998). Here, we evaluated link strength distri- 
butions based on the following two definitions: (i) the proportion of the prey in the predator's 
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Figure 4: The distribution of link strengths l^j^ > for the Lotka-Voherra model without 
adjusted scores (Fig. A) with adjusted scores (Fig. B) and the ratio-dependent functional 
response (Fig. C). The y— axes on all graphs are scaled logarithmically. The insets show the 
data with both axes scaled logarithmically, for Figs. B and C these include power-law fits 
which had exponents -0.74 and -0.29 respectively. 



diet after the population dynamics had equilibrated, 

4]' = ^. (8) 

and (ii) the per-capita interaction strength l\j . These are sometimes known as the elements 
of the "community matrix" (Laska & Wootton, 1998), although this phrase is also used to 
refer to the unnormalised matrix of partial derivatives (May 1973). It quantifies the strength 
of all direct interactions, both predator-prey interactions and interference competition. It is 
defined as: 

Definition (i) typically leads to a U-shaped link-strength distribution, with a peak at weak 
links and another one near the maximum link strength 1, indicating that many predators 
have one main prey. (Note that the second peak does not occur if the link strength is defined 
as the biomass passing through a link, since in this case the link strength is not normalised 
to X]j hj = !)• Definition (ii) typically leads to broad distributions with or without a peak at 
zero. 

Fig. m shows the distribution of link strengths Zfj^ for the two Lotka-Volterra models 
(without and with adjusted scores), and for the model with the ratio-dependent functional 
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Figure 5: The distribution of link strengths Z-^ > for the Lotka-Volterra models without 
adjusted scores (Fig. A) and with adjusted scores (Fig. B). 



response, averaged over many evolutionary time steps and several different model webs. For 
small link strength, the distribution for the ratio-dependent functional response resembles a 
power law with an exponent around —0.74, while the Beddington model with adjusted scores 
shows an exponent around —0.9, and the Lotka-Volterra model with adjusted scores shows an 
exponent around —0.29 (if a power law fit should be attempted at all). The Lotka-Volterra 
model without adjusted scores does not have several trophic layers, and its link strength 
distribution appears to decay exponentially at small values. For the Beddington functional 
response without adjusted scores, the decay is much steeper than with adjusted scores, but not 
exponential. All the above mentioned models are skewed towards small values except for the 
Lotka-Volterra model with adjusted scores where the approach to the origin is quite flat. For 
this model however only a small proportion, 3%, of links have l^j^ > 0. We conclude that for 
the models capable of generating large ecosystems only a small fraction of links were realised 
or llj^ had a large weight at small values. We can summarise this by calculating the fraction of 

links with zjj'' < 0.01 including zero, which gives 91% and 97% for the Lotka-Volterra models 
(without and with adjusted scores) and 97% for the ratio-dependent model. 

We can investigate these ideas further by examining the link strength according to the 
definition (ii); this can be either positive or negative. In the case of the Lotka-Volterra models, 

(2) 

the only direct inter-specific interactions are between predators and prey with llj = Xaij and 
/j^"* = —aij if i predates j. Thus in Fig. El we only need to show the positive half of the 

(2) 

llj distribution. From this we see that, whereas the Lotka-Volterra model without adjusted 
scores has a distribution of values that is skewed towards zero, the model with adjusted 
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Figure 6: The distribution of non-zero —aij for the Lotka-Volterra models with direct inter- 
specific competition. The results for the model without adjusted scores is shown in Fig. A 
and with adjusted scores in Fig. B. The fraction of inter-specific interaction with non-zero 
Oij was 100% for the former model and 23% for the latter. Note that the non-zero aij are 
clustered near to the smallest allowed magnitude of 2.5 in Fig. A and near to the largest 
allowed magnitude of 5 in Fig. B. 





scores has a maximum at an intermediate value of . However, as was the case for , we 
find a much higher proportion of non-zero links in the model without adjusted interactions 
(51% as compared to 3%). 

If we introduce direct inter-specific competition into the Lotka-Volterra models as in Sec- 
tion 3 by setting aij = cpij{l + qij)/'^ in Eq. (jT)), where pij is equal to one if i and j share at 
least one prey and is zero otherwise and where qij is the fraction of discrete features shared 
by the two species, we find that, even if adjusted interactions are used, only small webs can 
be evolved. The reason for this seems to be that the resulting direct competitive interactions 
which are allowed by the adjusted interactions are both strong and non-zero for a large frac- 
tion of species pairs in the web. This is illustrated by Fig. IHl where the distributions of the 
quantity —aij when c = 5.0 are plotted for the model with and without adjusted scores. We 
plot —aij because it is the contribution of direct competition to l\j . Note that because qij is 
discrete so is aij and that the range of non-zero aij is —5 < —aij < —2.5. 

The model with the ratio-dependent functional response, Eq. also includes direct inter- 
specific competition, as a result of the implicit sharing of prey between predators. However 
this does not lead to the strong competitive interactions seen in Fig. El probably because 
this competition occurs through the predator-prey terms, which as was shown in Fig. |3J are 
themselves very diffuse. This can be seen in Fig. where the distribution of if^ values for 
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Figure 7: The distribution of 1]^ for competitive interactions where the shared prey comprised 
greater than 1% of both competitors diet estimated for the ratio-dependent modeL This 
condition was satisfied by 5% of all species pairs. 



all competitive interactions and where the shared prey comprised greater than 1% of both 
competitors diet, is shown. This distribution, although bimodal, is heavily weighted towards 
the origin. This then explains why the ratio-dependent model with inter-specific competition 
is, in contrast to the Lotka-Volterra models discussed immediately above, capable of evolving 
large complex ecosystems. 

This finding of ours is complemented by a recent paper by Kondoh (2003), who investigates 
randomly linked food web models and the cascade model based on Lotka-Volterra dynamics 
with adjustable foraging efforts. He finds that incorporating adaptive foragers turns a negative 
relationship between stability and complexity, defined in terms of species number and potential 
connectance, into a positive one. This result had already been found, albeit over a more 
limited range of parameters, by Pelletier (2000). These results are strictly different to ours in 
that the stability criterion, community persistence defined as fraction of species surviving for 
some long time in a stochastic environment, is different. However the mechanism that allows 
the generation of large webs with high potential connectance, the reduction in the number of 
non-zero realised links and the reduction in strength of those links that are realised due to the 
foraging dynamics, is very similar to that proposed above. In fact this mechanism is arguably 
just May's original hypothesis, that strong links or a large fraction of links will destabilise 
ecosystems, placed in an adaptive context (May, 1972). 
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5 Discussion 



The form that a reahstic functional response might take has been the subject of a large number 
of papers in the literature but few, if any, of these suggestions were implemented in a model 
of multispecies communities in order to test their effectiveness. In this paper we have used an 
existing model to investigate the effect that different choices for the functional response have 
on food web structure. We studied models with a wide range of functional responses: Lotka- 
Volterra with and without direct inter-specific competition, HoUing and Beddington forms, 
all of these but with a mechanism to limit the number of species a predator actually feeds on, 
and ratio-dependent functional responses. We found that in the first two cases complex webs 
could not be built up, but in the last two cases they often could. More specifically, unless 
an (artificial) mechanism was introduced which restricted predator choice to prey which they 
were best suited to exploit or alternatively a ratio- dependent functional response was used, 
then stable webs consisting of more than one trophic level could not be built up. 

Our second major conclusion was that, given we chose a functional response which gave 
rise to complex, stable webs as described above, the link strength distributions were skewed to- 
wards zero. For the Beddington functional response with adjusted scores or ratio-dependence, 
the link strength distribution followed a power-law, but this was much less clear in the case 
of the Lotka-Volterra model, where only a small fraction of the links were realised i.e. a large 
fraction of the links were zero. Some insight into the mechanisms involved in these effects was 
also gained by looking at the strength of a direct inter-specific competition introduced into 
the Lotka-Volterra equations. Without adjusted scores, interactions were present, but not 
particularly strong. However with adjusted scores, about three-quarters of the interactions 
(for the parameter values we used) were zero, but the remaining quarter were strong. In this 
case we found that, even with adjusted scores, large complex webs could not be grown. We 
concluded that the existence of strong competitive interactions might also destabilise food 
webs. 

In summary, we found that the type of functional response used in the population dynamics 
of multispecies communities has to be chosen carefully if a large complex community is to be 
sustained. If this is achieved, then a large proportion of weak links arises naturally from the 
evolution of the food webs. These conclusions were arrived at within the context of a class of 
evolutionary food web models. It would be interesting to investigate if these results could be 
obtained from other starting points or with different model assumptions. 
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